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ABSTRACT 


J  An  experiment  was  performed  to  study  the  effects  of  an 
underwater  explosion  on  a  submerged  test  panel.  This  is  a 
part  of  an  on-going  research  project  sponsored  by  the  Defense 
Nuclear  Agency  (DNA)  to  help  predict  the  dynamic  responses 
of  naval  structural  members  due  to  underwater  shocks.  An 
important  element  of  this  thesis  was  to  validate  a  finite 
element/finite  central  difference  computer  code  developed  to 
forecast  shell  responses.  Emphasis  was  placed  upon  attaining 
stiffener  tripping  and  collecting  invaluable  dynamic  flat 
plate  tripping  data.— .  <  ■< ''  *■  ■  1 
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I .  INTRODUCTION 


A.  OBJECTIVES 

Current  building  specifications  for  U.S.  Navy  ships  and 
submarines  require  shock  to  be  met  on  shipboard  equipment 
and  that  underwater  shock  testing  be  performed  on  at  least 
one  ship  in  each  class  [Ref.  1] .  Underwater  explosion 
(UNDEX)  research  is  conducted  at  NPS  because  of  the  impor¬ 
tance  of  shock  testing  to  the  Navy.  The  objectives  of  this 
thesis  are  four-fold.  First,  to  investigate  the  large  de¬ 
flection  elasto-plastic  transient  response  characteristics  of 
submerged  stiffened  flat  plates,  subjected  to  underwater  shock. 
Secondly,  to  perform  underwater  explosion  testing  using  a 
stiffened  flat  plate  and  to  validate  the  predicted  results 
with  those  of  the  experiment.  Thirdly,  to  investigate  stiffener 
tripping  and  its  effect  to  the  gross  shell  responses.  And 
finally,  to  acquire  the  knowledge  and  technical  expertise 
necessary  to  conduct  these  types  of  experiments. 

What  is  stiffener  tripping  and  why  is  it  significant  to 
naval  construction?  Stiffener  tripping  is  defined  as  the 
lateral-torsional  instability  of  a  stiffener  as  it  becomes 
suddenly  unstable  and  fails  under  load  [Ref.  2].  Tripping 
is  identified  by  a  twisting  of  the  stiffener  about  its  line 
of  attachment  to  a  structure.  The  purpose  of  stiffeners  is 
to  add  strength  to  ship  hulls.  Removal  of  the  stiffener 


weakens  the  structure.  A  report  published  in  1979  by  John 

Adamchak  stresses  the  importance  of  stiffener  tripping 

effects  on  shipboard  hulls  and  grillages  by  stating  that 

previous  underwater  tests 

clearly  demonstrate  the  significance  of  lateral- 
torsional  instability  (tripping)  as  a  primary  ductile 
failure  mode  for  a  ship's  structure.  [Ref.  2] 

Once  a  stiffened  structure  has  experienced  tripping,  its 

ability  to  continue  to  resist  deformation  under  loading  is 

drastically  reduced. 

B .  BACKGROUND 

This  is  the  third  NPS  thesis  in  the  field  of  UNDEX  and 
their  dynamic  effects  on  submerged  plates.  This  research  is 
part  of  a  major  project  coordinated  and  sponsored  by  the 
Defense  Nuclear  Agency.  During  the  early  part  of  19  84  ,  Lt . 
Tom  Rentz  conducted  the  first  NPS  thesis  research  on  UNDEX 
responses  of  flat  plates.  A  plate  positioned  on  top  of  an 
air-backed  test  platform  was  chosen  to  best  simulate  the 
conditions  present  in  ship  hulls. 

Fabrication  and  exact  specifications  of  the  test  plat¬ 
form  are  covered  in  detail  in  [Ref.  3],  with  Figure  1.1 
showing  a  sketch  of  this  structure.  Flat  plates  are  manu¬ 
factured  out  of  6061-T6  aluminum  because  it  is  accessible, 
easy  to  machine  and  this  material  exhibits  strain  rate  inde¬ 
pendent  properties.  The  first  plate  tested  is  shown  in 
Figure  1.1  and  consisted  of  two  parallel  rectangular  bar 
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1*5**  (TYP) 


stiffeners  with  a  thickness  of  0.1875  inches  and  height  of 
1.0  inches  running  transversely  on  the  plate. 

In  February  1984,  an  explosive  shot  was  performed  on  the 
first  stiffened  plate.  Last-minute  changes  made  it  necessary 
for  sixteen  half-pound  rectangular  blocks  of  TNT  to  be  sub¬ 
stituted  for  a  single  eight-pound  charge.  Test  results 
revealed  large  deformations  and  tears  at  stiffener  ends 
which  prevented  stiffener  tripping. 

Lt.  Nelson  King  was  the  next  to  pursue  the  area  of  UNDEX 
at  NPS .  During  the  summer  months  of  1984,  he  performed  two 
explosions,  using  different  plate  designs  for  each.  His 
first  plate  was  similar  to  Figure  1.1  except  that  the 
stiffeners  were  reconfigured  to  a  tee,  see  Figure  1.2.  It 
was  believed  by  exposing  more  surface  area  to  the  shock 


A 


Figure  1.2  King's  Tee-shaped  Stiffener 


wave,  a  larger  compression  zone  would  be  developed,  encourag¬ 
ing  tripping.  A  single  nine  and  a  half-pound  test  charge 
was  exploded.  Test  results  were  similar  to  the  first  explo¬ 
sion.  A  move  towards  a  simpler  test  panel  was  taken. 

The  test  panel  was  redesigned  to  a  single  bar  stiffener 
with  increased  thickness  to  0.25  inches  and  height  to  1.25 
inches  running  longitudinally.  Figure  1.3  shows  the  geometry 
of  this  new  panel.  A  half-pound  charge  of  TNT  was  exploded 
over  the  test  plate.  Recorded  sensor  results  were  invalid 
for  this  experiment  but  micrometer  measurements  taken  later 
indicated  that  perhaps  a  slight  amount  of  stiffener  bowing 
was  taking  place  and  further  testing  with  this  test  geometry 
was  warranted  [Ref.  4]  . 

The  test  panel  investigated  for  this  thesis  is  identical 
to  the  single  longitudinal  stiffener  model.  A  charge  weight 
of  eight  pounds  TNT,  a  depth  of  four  feet  and  a  plate  stand¬ 
off  distance  of  nine  feet  was  selected.  The  experiment  was 
conducted  at  West  Coast  Shock  Facility  (WCSF)  at  Hunter's 
Point  Shipyard,  California. 


Figure  1.3  King's  Plate  and  Current  Model 


II.  UNDERWATER  EXPLOSION  THEORY 

Before  continuing  with  the  experimental  aspect  of  this 
research,  an  overview  of  UNDEX  theory  is  provided.  Explosions 
are  a  very  complex  field  of  study.  The  destructive  effect 
of  an  explosion  is  dependent  upon  the  composition  and  geometry 
of  the  explosive  element,  the  media  or  surrounding  environ¬ 
ment,  the  stand-off  distance  to  the  target  and  the  physical 
nature  of  the  target.  Simply  stated,  an  explosion  is  a  chemi¬ 
cal  reaction  which  releases  energy;  both  thermal  and  mechanical 
[Ref.  5] . 


Figure  2.1  Shock  Wave  Pressure  Profile 


A.  THE  SHOCK  WAVE 


An  immediate  result  of  underwater  detonation  is  a  high 
gas  pressure  build-up.  The  gases  interact  with  the  surround¬ 
ing  water  causing  compression  and  a  radial  outward  flow. 

This  forms  the  shock  wave  which  comprises  about  half  of  the 
explosion's  energy.  It  is  the  effects  of  the  shock  wave  on 
the  flat  plate  that  is  the  primary  area  of  interest  for  this 
thesis  research.  The  exact  character  of  the  shock  wave  near 
the  charge  is  not  well  known,  thus,  focus  will  be  on  the 
characteristics  of  the  shock  wave  at  some  distance  away. 
Experience  has  shown  that  the  shock  wave  travels  at  the  speed 
of  sound  in  water  (approximately  5000  feet  per  second)  [Ref. 

6] .  Formulas  have  been  developed  for  distances  between  10 
and  100  times  the  charge's  radius  to  describe  such  features 
as:  peak  pressure,  pressure  as  a  function  of  time,  and  time 
decay  constant.  Ideally,  the  shock  wave,  which  delivers  a 
pressure  incident,  has  a  profile  similar  to  Figure  2.1.  There 
is  a  sharp  peak  pressure  spike,  followed  by  an  exponential 
decay,  with  about  one  third  the  original  pressure  occurring 
after  only  a  single  time  decay  constant. 

A  few  equations  useful  in  the  description  of  the  shock 
wave  pressure  profile  are  listed  below.  Pressure  is  des¬ 
cribed  by  equation  2.1. 

(t-t^ 

"  5 

P(t)  =  P max  e  ;  given  t  >  t^  psi  (2.1) 


0  is  the  explosive  decay  constant  and  P 

c  1  max 

given  by  equations  2.2  and  2.3. 
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The  time  of  shock  wave  arrival  at  the  target,  T,  is  calcu¬ 
lated  from  equation  2.4. 


T 


R 

C 


sec 


(2.4) 


where : 

W  is  explosive  weight  in  pounds 

R  is  stand-off  in  feet  from  charge  to  target 

t-t^  is  elapsed  time  after  shock's  arrival 

A-^,  A2,  and  K2  are  constants  which  vary  depending 
on  the  type  of  explosive  (see  Table  I) 

C  is  the  speed  of  sound  in  water 


The  experimental  values  expected  for  an  eight-pound  charge  of 
TNT  with  a  nine  foot  target  stand-off  are  peak  pressure  of 
3814  psi,  decay  constant  of  0.153  msec  and  shock  arrival  time 


of  1.8  msec. 


TABLE  I 


Explosive 

Constants 

Explosive 

TNT 

HBX-1 

Pentolite 

K1 

22,505 

22,347.6 

24,589 

A1 

1.18 

1.144 

1.194 

K2 

.058 

.056 

.052 

a2 

-.185 

-.247 

-.257 

K3 

12.67 

14.14 

12.88 

K4 

4.268 

4.761 

4.339 

In  order  to  obtain  clean  data  the  secondary  effects  of 
UNDEX  had  to  be  eliminated  or  reduced.  UNDEX  secondary 
effects  to  a  target  include  gas  bubble,  bulk  cavitation,  and 
surface  and  bottom  reflection  effects.  A  description  of 
each  of  these  phenomena  is  provided. 

B.  THE  GAS  BUBBLE 

The  gas  bubble  is  an  UNDEX  phenomenon  which  is  closely 
related  to  the  shock  wave.  At  detonation,  while  the  shock 
wave  flows  radially  outward,  the  highly  compressed  gases 
form  a  bubble  which  expands  in  size  with  an  accompanying 
decrease  of  bubble  gas  pressure.  The  internal  gas  pressure 
gradually  decays  as  the  bubble  expands  until  hydrostatic 
pressure  finally  takes  over  forcing  the  bubble  to  contract. 
Bubble  contraction  is  slow  at  first  but  progresses  rapidly 
until  a  minimum  bubble  radius  is  attained.  The  bubble  even¬ 
tually  collapses  producing  a  pressure  pulse  which  stimulates 


radial  flow  outward  forming  a  new  bubble  that  initiates  the 
cycle  all  over  again.  Pressure  pulses  continue  to  be  pro¬ 
duced  until  all  of  the  bubble  energy  is  emitted.  These 
pulses  make  up  part  of  the  reloading  nature  of  an  underwater 
explosion.  Although  the  bubble  pressure  pulses  are  signifi¬ 
cantly  weaker  in  magnitude  than  the  incident  shock  wave  (on 
the  order  of  10-15  percent  of  the  first  bubble  pulse) ,  they 
still  represent  important  dynamic  loads  to  a  target  [Ref.  7], 
An  interesting  aspect  of  gas  bubble  theory  is  the  migra¬ 
tion  of  the  bubble  to  the  surface.  The  migration  is  upward 
due  to  the  influence  of  gravity,  with  the  most  rapid  changes 
noticed  when  the  bubble's  volume  is  at  a  minimum.  Figure  2.2 
depicts  the  relationship  between  bubble  expansion,  pressure 
pulse  emissions  and  bubble  migration.  This  figure  is  adapted 
from  [Ref.  7] .  Experiments  have  shown  that  the  energy  asso¬ 
ciated  with  the  gas  bubble  is  about  41  percent  of  the  total 
energy  available  during  an  explosion  [Ref.  1] .  The  following 
empirical  formulas  are  presented  to  give  some  appreciation 
for  the  gas  bubble  and  its  behavior.  The  maximum  bubble 
radius  and  the  period  of  the  first  pressure  pulse  are  given 
by  equations  2.5  and  2.6. 


3  (D+33)1/3 


(2.5) 


T  =  K 
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Figure  2.2  Bubble  Historical  Representation 


whe  re : 


is  maximum  bubble  radius 

max 

K,  and  K.  are  constants  from  Table  I 
3  4 

D  is  the  submerged  charge  depth  in  feet 

Both  of  these  equations  employ  the  fact  that  the  pressure  at 
the  surface  of  the  water  is  one  atmosphere  or  33  feet  of  sea 
water. 

• 

C.  BULK  CAVITATION  ZONE 

Bulk  cavitation  is  an  area  where  air  bubbles  or  pockets 
have  formed  which  are  a  result  of  the  reflection  of  the 
compressive  shock  wave  off  the  water's  surface.  As  the 
shock  wave  reflects  from  a  free  surface  it  changes  in  sign 
to  a  tension  wave.  The  water  mass  eventually  pulls  apart 
beneath  the  free  surface  to  prevent  any  further  tension  and 
voids  are  created  [Ref.  6] .  The  parameters  which  affect  the 
area  of  the  bulk  cavitation  zone  are  charge  type,  weight  and 
depth . 

D.  FREE  SURFACE  AND  BOTTOM  REFLECTIONS 

Free  surface  and  bottom  reflections  are  termed  'late¬ 
time'  effects.  The  shock  wave  moves  radially  outward  from 
the  charge,  very  much  like  the  resultant  water  motion  from  a 
pebble  thrown  into  a  still  pond.  The  incident  shock  wave 
proceeds  directly  to  the  target,  while  the  bottom  and  free 
surface  serve  as  boundaries  or  obstructions.  Surface  cut-off 


occurs  at  the  target  when  the  reflected  free  surface  negative 
pressure  wave  strikes.  To  help  visualize  what  is  happening 
with  surface  reflection,  an  image  source  is  used.  Figure  2.3 
shows  the  assumed  image  source,  charge  and  point  of  interest/ 
target.  Figure  2.4  depicts  the  type  of  pressure  profile 
experienced  by  the  target-.  Using  the  geometry  of  Figure  2.3, 
surface  cut-off  time  and  pressure  can  be  computed  by  equations 
2.7  and  2.8.  Both  of  these  figures  are  taken  from  [Ref.  1]. 
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Figure  2.4  Incident  and  Surface  Reflected  Pressure  Waves 
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The  influences  of  the  sea  bed  or  bottom  are  not  as  well 


understood  as  those  of  the  surface.  However,  the  bottom  is 
considered  to  be  a  perfectly  rigid  boundary  reflecting  waves 


E.  EXPERIMENTAL  CONDITIONS  AND  PRE-SHOT  CALCULATIONS 

As  mentioned  earlier  one  of  the  objectives  of  this  thesis 
was  to  obtain  clean  data,  looking  primarily  at  the  effects 
of  the  incident  shock  wave.  By  careful  attention  to  test 
platform  positioning  in  the  water,  most  of  the  secondary 
effects  could  be  eliminated.  A  charge  depth  of  four  feet  was 
selected  which  allowed  for  adequate  venting  of  the  gas 
bubble  to  the  atmosphere  prior  to  hitting  the  test  plate. 
Additionally,  experience  has  shown  that  a  satisfactory  vent¬ 
ing  depth  is  one  which  is  between  50  and  80  percent  of  maximum 
gas  bubble  radius  [Ref.  1].  Using  equation  2.5,  a  maximum 
bubble  radius  of  7.6  feet  was  computed.  This  represents  a 
depth  to  radius  ratio  of  52.6%. 

To  eliminate  bulk  cavitation,  the  cavitation  region  had 
to  be  first  determined.  The  equations  necessary  to  calculate 
this  region  are  extensive  and  tedious,  thus  a  computer  program 
was  used  to  graph  the  cavitated  region,  see  Figure  2.6.  To 
avoid  bulk  cavitation  and  in  keeping  with  past  NPS  UNDEX 
experiments,  a  nine  feet  stand-off  distance  was  used. 

Bottom  and  surface  reflections  are  late  time  effects. 

By  positioning  the  test  structure  directly  beneath  the  charge, 
there  is  a  time  window  of  approximately  2.5  msec  where 
reflections  can  be  ignored.  The  experiments'  sensors  were 
sampled  during  this  time  frame.  The  depth  of  the  water  at 
the  test  site  is  approximately  thirty  feet. 
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Figure  2.6  Bulk  Cavitation  Zone 
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III.  THE  COMPUTER  CODE 


The  computer  code  selected  for  use  in  this  thesis  is 
EPSA,  Elasto-Plastic  Shell  Analysis,  developed  by  Weidlinger 
Associates.  This  is  the  second  attempt  at  NPS  to  model  a 
flat  plate  based  on  the  EPSA  code.  Lt .  Rentz  was  the  first 
to  utilize  EPSA,  see  [Ref.  3], 

EPSA  is  a  finite  element,  finite  central  difference  com¬ 
puter  code.  For  calculation  purposes  EPSA  assumes  a  symmetric 
loading  pattern  and  time  zero  commences  when  the  spherical 
shock  wave  first  encounters  the  test  plate.  DNA/NAVS EA  funded 
its  development  for  the  purpose  of  analyzing  dynamic  loading 
effects  on  submerged  stiffened  shells.  The  theory  incorporated 
in  EPSA  is  explained  in  the  EPSA  User's  Manual  [Ref.  8], 

The  purpose  of  this  chapter  is  to  describe  how  EPSA  was 
used  in  the  pre-shot  calculations  for  this  experiment.  In 
Chapter  V,  EPSA  will  be  reviewed  emphasizing  its  post-experiment 
capabilities.  EPSA  is  operated  on  a  Vax-11/780  computer 
system  located  in  Spanagel  Hall. 

Symmetry  conditions  were  incorporated  to  model  the  aluminum 
test  plate  for  EPSA,  see  Figure  3.1.  Requiring  only  half 
of  the  plate  to  be  discretized  saved  on  computer  time  and 
storage  space.  EPSA  computes  quantities  from  nodal  points 
and  elements.  Plate  displacements  and  pressures  are  nodal 
point  values,  whereas  stresses  and  strains  are  elemental. 


Figure  3.1  Half  Plate  EPSA  Discretization  Model 
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Selecting  a-  discretization  of  fifteen  columns  by  twenty  five 
rows  yielded  rectangular  block  sizes  of  0.6  by  0.48  inches. 

The  total  number  of  elements  and  nodes  are  375  and  416, 
respectively.  Sides  one,  three,  and  four  have  fixed  boundary 
conditions  and  side  two  represents  a  symmetry  edge. 

A  computer  data  deck  is  necessary  to  run  the  EPSA  code. 

A  sample  data  deck  created  for  this  thesis  is  reproduced  in 
Figure  3.2.  The  User's  Manual  describes  in  detail  the  infor¬ 
mation  included  in  this  input  deck.  It  is  through  the  input 
deck  that  the  plate  physical  characteristics  and  the  under¬ 
water  experimental  conditions  are  passed  to  the  computer  code. 

To  facilitate  in  the  viewing  of  EPSA  output  files,  Francois 
Daube,  a  thesis  student,  wrote  an  interface  program  with 
PATRAN-G.  PATRAN  provides  a  color  graphics  capability  for 
the  display  of  von  Mises'  stresses  and  plate  displacements 
at  requested  time  increments.  Daube' s  work  is  presented  in 
[Ref.  10]  and  the  EPSA-PATRAN  interface  is  covered  in  [Ref.  3]. 

A  color  chronological  history  of  von  Mises'  stresses  and 
displacements  for  an  eight-pound  TNT  charge  was  created. 

By  looking  at  PATRAN  and  EPSA, pre-shot  outputs  strain  gage 
locations  were  selected.  See  Figure  3.3.  Emphasis  was 
placed  upon  positioning  the  gages  in  an  X,Y  orientation  and 
maintaining  symmetry  in  order  to  compare  experimental  load¬ 
ing  patterns  to  EPSA  predictions. 
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Figure  3.2  A  Sample  EPSA  Computer  Deck 
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IV.  UNDERWATER  SHOCK  TESTING 


A.  INSTRUMENTATION 

The  UNDEX  experiment  was  conducted  at  the  West  Coast 
Shock  Facility  (WCSF) .  This  is  a  Navy  sponsored  activity 
accessible  to  NPS ,  located  about  120  miles  to  the  north, 
in  the  South  Bay  of  San  Francisco. 

Thirteen  strain  gages,  two  accelerometers  and  two  pressure 
transducers  were  installed  on  the  test  structure.  The  loca¬ 
tion  of  the  strain  gages  has  already  been  shown  on  Figure  3.3. 
The  accelerometers  were  mounted  alongside  each  other  at  the 
bottom  of  the  test  structure.  One  of  the  pressure  transducers 
was  placed  on  a  nine  foot  arc  from  the  charge  approximately 
five  feet  away  from  the  test  platform  to  record  incident 
pressure  data.  This  was  secured  in  place  by  light  manila 
guidelines  and  anchored  vertically  by  a  weight.  The  other 
was  situated  on  top  of  the  plate  coming  out  of  one  of  the 
corners  at  a  45°  angle  towards  the  plate's  center.  This 
transducer  served  as  the  basis  for  the  plate's  total  pressure 
measurements.  The  sensor  specifications  are  listed  in  Table 
II. 

Two  Honeywell  M101  Wideband  II  (direct  record)  tape  units 
operating  at  120  inches-per-second  recorded  all  data  channels. 
On-site  verification  of  data  was  accomplished  with  a  Honeywell 
1508b  Visicorder.  Post-processing  of  the  data  was  through 


TABLE  II 


Instrumentation  Specifications 


EQUIPMENT 
Strain  gages 
Pressure  gages 


TYPE 

CEA-350  ohms 
.25"  Tourmaline 


MAX  RANGE 

50,000  microstrains 

10  kpsi,  97%  Response 
Ratio 


Accelerometers  PCB  Piezoelectric  2500  and  100,000  G's 
Amplifiers  Ektron  563FJ  - 


an  HP-5451C  Fourier  Analyzer  located  at  NPS.  A  schematic  of 
the  instrumentation  is  shown  in  Figure  4.1.  To  assist  in 
data  reviewing,  an  ingenious  triggering  mechanism  was  devised 
by  electronic  technician  Tom  Christian.  This  triggering 
device  consisted  of  a  small  DC  circuit  wrapped  around  the 
charge.  At  time  of  detonation  the  wire  breaks  registering 
a  zero  voltage  reading  on  the  recorder. 

Instrumentation  wires  ran  down  from  the  connection  box 
on  the  test  structure,  along  the  bay's  floor  to  the  WCSF's 
instrumentation  shack.  Figure  4.2  is  a  photo  of  the  Honey¬ 
well  recorders  and  their  associated  amplifiers  in  the 
electronic's  shack.  It  should  be  noted  that  a  redesigned 
connection  box  was  constructed  and  implemented  for  this  shot. 
On  the  previous  UNDEX  experiments  a  brass  electrical  box  was 
installed,  which  did  not  afford  adequate  protection  to  the 
connections.  To  remedy  this  weakness,  a  new  box  was  fabri¬ 
cated  out  of  1.25  inch,  A36  structural  steel  plate. 


B.  SHOT  DAY 


Thursday,  6  December  1984,  was  set  as  the  shot  date.  A 
drawing  of  the  test  rig  is  provided  by  Figure  4.3.  The 


Figure  4.3  Test  Rig  Model 

weather  was  ideal  for  an  UNDEX  experiment — clear  and  sunny 
skies,  calm  seas  and  no  wind.  After  re-examining  the  equip¬ 
ment  and  performing  calibration  checks  on  the  sensors,  the 
riggers  took  over  preparing  the  test  structure  for  the  bay. 
Figures  4.4,  4.5,  and  4.6  show  the  shot  sequence.  Immediately 
after  the  explosion,  everyone  sensed  something  had  gone  wrong 
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with  the  shot  as  the  water  plume  was  smaller  than  expected, 
and  the  explosive  sound  was  muffled.  Tom  Christian,  stationed 
in  the  instrumentation  shack,  yelled  out  moments  later  that 
he  had  picked  up  some  very  weak  signals  on  all  channels,  but 
was  unsure  of  their  meaning.  Figure  4.7  reveals  the  test 
structure  as  it  is  being  hauled  out  of  the  water  following 
the  explosion.  No  significant  external  damage  was  apparent 
and  it  was  noticed  that  the  free  field  pressure  transducer 
was  intact  and  still  properly  oriented.  It  was  not  until  the 
test  structure  was  repositioned  on  dock,  that  a  large  amount 
of  unexploded  TNT  was  discovered  on  top  of  the  plate,  see 
Figure  4.3.  The  appearance  of  the  TNT  on  the  plate  was  con¬ 
clusive  that  the  full  eight-pounds  of  explosive  did  not  go 
off.  The  unexploded  residue  was  collected  and  weighed  approxi¬ 
mately  one  and  a  half  pounds.  During  post-analysis,  by 
examining  the  pressure  graphs  from  the  Fourier  Analyzer,  it 
was  determined  that  approximately  0.16  pounds  of  TNT  actually 
exploded.  This  was  obtained  by  inputting  the  recorded  value 
of  820  psi  into  equation  2.3  and  solving  for  charge  weight. 

The  plate  surface  was  carefully  inspected  with  no  visible 
damages  noted.  Deflection  measurements  were  not  taken  as  a 
micrometer  was  not  available.  The  electrical  connection  box 
remained  water-tight  and  all  but  one  of  the  sensors  were 
still  operational.  After  a  short  deliberation  and  determin¬ 
ing  the  availability  of  another  eight-pound  charge,  it  was 
decided  to  fire  another  shot  later  that  same  day.  WCSF 


personnel  had  stated  that  this  was  the  first  incomplete 
explosion  of  this  type  they  had  experienced.  Incidentally, 
it  should  be  pointed  out  that  firing  a  small  primer  charge 
before  a  larger  explosion  is  not  an  uncommon  practice. 
Experience  has  shown  that  a  primer  tends  to  aid  in  gage 
calibration  and  loosens  up  residual  plate  stresses. 

Approximately  three  hours  later,  another  explosion 
occurred.  This  time  two  blasting  caps  were  used  to  ensure 
full  detonation  of  the  TNT.  It  was  apparent  from  the  noise 
and  water  plume  that  the  second  shot  was  significantly 
larger  them  the  early  one.  It  was  later  determined  from 
pressure  data  graphs  that  the  peak  pressure  pulse  experienced 
by  the  plate  for  this  second  shot  was  4217  psi,  which  is 
equivalent  to  a  9.84  pound  TNT  charge.  Figure  4.9  shows  the 
water  plume  produced  by  the  second  shot.  One  of  the  pneumatic 
bladders  ruptured  from  the  force  of  the  explosion,  as  shown 
in  Figure  4.10.  Another  post  shot  view.  Figure  4.11,  shows 
the  damage  experienced  by  the  test  rig.  Notice  how  two  of 
the  test  platform  support  cables  have  been  torn  free. 

A  view  of  the  underside  of  the  test  structure  and  cable 
runs  going  into  the  electrical  connection  box  is  shown  in 
Figure  4.12.  In  the  past  the  electrical  cables  have  been 
pulled  loose  from  the  box,  exposing  the  connections  to  the 
environment.  Following  this  second  explosion,  Tom  Christian 
joyfully  announced  that  data  had  been  collected. 
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V.  RESULTS  AND  RECOMMENDATIONS 

A.  PHYSICAL  RESULTS 

As  mentioned  in  Chapter  IV,  the  first  shot  did  not  pro¬ 
duce  any  noticeable  physical  results.  However,  the  second 
shot  did  achieve  large  deformations  as  evidenced  by  Figures 
5.1  and  5.2.  Both  stiffener's  ends  sheared  and  there  are 
tears  along  three  of  the  sides.  The  exact  tear  lengths  are 
recorded  on  Figure  5.3.  Not  visible  in  any  of  the  photos  is 
the  fact  that  several  strain  gages  became  detached  from  the 
plate  surface  during  the  testing.  Additionally,  a  few  of  the 
strain  gage  lead  wires  were  severed.  This  severing  was 
probably  due  to  a  guide  wire  shackle  that  parted  from  the 
ruptured  flotation  device  following  the  shot. 

Dial  indicator  deflection  readings  were  taken  while  the 
aluminum  plate  was  still  clamped  to  the  test  structure,  and 
are  presented  in  Figure  5.3.  An  EPSA  run  was  made  using  a 
9.84  pound  TNT  charge  to  determine  predicted  deflections. 

EPSA  predicted  and  actual  displacements  are  compared  on 
Figure  5.3.  There  are  large  discrepancies  noticed  in  these 
comparisons,  especially  near  locations  where  plate  tears 
exist.  EPSA  does  not  take  into  account  possible  secondary 
effects  which  may  have  loaded  the  plate,  most  notably  the 
loading  generated  by  water  particles  striking  the  plate  due 
to  the  downward  rigid  body  action  of  the  test  structure. 


gure  5.1  Shearing  at  Both  Stiffener  Ends 


Figure  5.2  Large  Plate  Tear  at  the  Side 
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EPSA  assumes  a  perfectly  symmetric  loading  pattern,  without 
any  stress  concentrators  and  does  not  incorporate  fracture 
mechanics . 


B.  DATA  INTERPRETATION 

The  majority  of  data  interpretation  was  conducted  at  NPS 
using  the  HP-5451C  Fourier  Analyzer.  Sensor  data  recorded 
on  FM  tape  was  transferred  from  the  Honeywell  recorder  to  a 
5451C  disk.  All  channel  data  was  recorded  at  120  inches-per 
second  (ips)  and  played  back  to  the  disk  at  1.875  ips.  A 
real  time  conversion  factor  of  64  is  necessary  to  change  from 
120  to  1.875  ips.  Time  zero  for  data  reduction  purposes  is 
the  exact  moment  of  detonation.  The  triggering  channel  was 
instrumental  and  served  as  a  reference  in  transferring  data 
from  the  recorder  to  the  disk.  All  seventeen  channels  are 
permanently  stored  both  on  tape  and  disk  at  NPS.  Each  data 
channel  consists  of  fifteen  storage  records,  with  each  record 
cracking  20.47  msec  at  1.875  ips  time  or  0.319  msec  real  time 

Table  III  gives  a  summary  of  data,  storage  locations  and 
conversion  scale  factors  used  for  shots  one  and  two.  Sample 
digitized  output  data  graphs  are  provided  in  Appendix  A. 

Time  was  not  available  to  write  a  computer  program  to  change 
the  scales  on  the  axes  to  reflect  proper  units  or  actual 
times.  However,  reference  values  are  indicated  on  each 
graph . 

As  mentioned  earlier,  the  digitized  incident  pressure 
data  from  shot  one  yielded  a  peak  pressure  of  820  psi, 
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TABLE  III  (CONTINUED) 


\ 
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equating  to  0.16  pounds  of  TNT  at  a  real  time  of  2.05  msec. 
This  value  was  compared  to  that  of  the  plate's  surface  total 
pressure  gage  of  585  psi  at  time  2.13  msec.  It  is  difficult 
to  say  why  the  peak  surface  pressure  is  less  than  the  peak 
incident  pressure  but  perhaps  the  time  difference  of  0.079 
msec  accounts  for  this  degree  of  decay.  This  statement  is 
predicated  on  the  assumption  that  both  pressure  transducers 
were  the  same  distance  away  from  the  charge.  Selecting  the 
value  of  0.16  pounds  of  TNT,  another  EPSA  computer  run  was 
made  to  determine  the  impact  of  this  weight.  Examining  the 
EPSA  output  indicates  the  possibility  that  some  permanent 
plate  deflections  may  have  resulted  from  the  first  shot. 
Maximum  center  plate  deflection  was  calculated  to  be  0.0147 
inches.  Using  the  aluminum  plate  thickness  of  0.25  inches, 
a  ratio  of  deflection  to  plate  thickness  of  4%  is  attained. 
This  is  considered  a  small  deflection  and  serves  as  a  basis 
to  ignore  the  impact  of  the  first  explosion.  Another  test 
performed  to  satisfy  the  hypothesis  that  the  first  shot  did 
not  have  an  overpowering  effect  on  the  test  plate  was  a 
comparison  of  peak  strains  to  yield  strain.  The  typical 
yield  strain  for  6061-T6  aluminum  is  4000  microstrains.  The 
EPSA  strains  produced  from  a  0.016  pound  TNT  load  are  all 
well  below  yield.  This  indicates  that  the  plate  did  not 
enter  its  plastic  region. 

Peak  strain  comparisons  are  provided  in  Figure  5.4.  The 
values  within  the  parentheses  are  EPSA  predicted  and  the 


Comparisons  of  Predicted  to  Actual  Strains 


ones  outside  represent  experimental.  It  was  disappointing- 
to  find  such  -large  discrepancies  between  gages.  The  measured 
peak  values  are  haphazard,  without  any  regard  for  consistency 
in  reading  higher  or  lower  than  expected.  Two  of  the  strain 
gages  gave  no  output  due  to  severed  lead  wires,  and  one  indi¬ 
cated  a  reading  beyond  gage  limits.  Accelerometer  magnitudes 
were  equally  useless.  Both  transducers  registered  values 
exceeding  their  maximum  range.  Sensor  symmetry  responses  were 
investigated  by  comparing  the  outputs  from  three  sets  of 
strain  gages  and  both  accelerometers.  Table  IV  gives  the 
highlights  of  these  comparisons.  As  alluded  to  earlier,  the 
peak  strain  values  are  inconsistent  but  loading  times  are 
very  close.  Poor  peak  value  comparisons  add  a  great  deal  of 
uncertainty  about  the  validity  of  the  strain  gage  measurements. 

TABLE  IV 

Symmetry  Sensor  Comparison 
SENSOR  SHOT  #1  SHOT  #2 

SG1&7  No  comparisons  Arrival  Times  within  0.7% 

due  to  bad  SGl  values  Peaks  off  by  factor  of  2.2 

SG3&6  Arrival  times  within  1.0%  Arrival  times  within  0.7% 

Peaks  off  by  factor  2.1  no  comparison,  bad  SG6  value 

SG10&13  not  recorded  Arrival  times  within  2.4% 

Peaks  off  by  factor  of  62 

X"  Arrival  times  within  0.3%  Arrival  times  within  2.7% 

Peaks  off  by  factor  1.16  Peaks  off  by  factor  1.4 

Notes:  Peaks  represent  strains  for  the  strain  gages 

and  G's  for  the  Accelerometers  (X") 


The  data  pressure  histories  from  shot  two,  enclosed  in 
Appendix  A,  seem  to  be  the  most  useful  data  obtained  from 


this  experiment.  The  peak  surface  total  pressure  sensor 
value  of  4.922  kpsi  is  larger  than  the  free  field  peak  of 
4.217  kpsi.  Theoretically,  plate  total  pressure  should  be 
twice  as  large  as  free  field;  but  this  is  assuming  an  ideal 
rigid  body  response.  An  aluminum  plate  does  not  reflect 
perfectly,  so  any  magnification  between  one  and  two  is 
practical . 

There  does  exist  a  time  difference  between  pressure 
transducers.  Prior  to  the  structure's  water  entry,  there  was 
exactly  a  nine  foot  charge  to  transducer  separation  distance. 
A  time  comparison  attempt  was  made  to  check  whether  these 
distances  remained  the  same  once  submerged.  The  results 
of  this  comparison  are  shown  in  Table  V.  Theoretically,  the 


TABLE  V 


Pressure  Sensor  Time  Comparisons 

SHOT  #1  SHOT  #2 


SENSOR 

ARRIVAL 

PEAK 

ARRIVAL 

PEAK 

P(t) 

2.116 

msec 

2.132 

msec 

2.390 

msec 

2.403 

msec 

P(f> 

1.950 

msec 

2.024 

msec 

2.185 

msec 

2.213 

msec 

Time  change 

0.166 

msec 

0.109 

msec 

0.203 

msec 

0.190 

msec 

Distance 

0.83 

ft 

0.54 

ft 

1.02 

ft 

0.95 

ft 

Note:  Distance  is  computed  assuming  the  speed  of  sound, 

5  ft/msec 


shock  wave  should  have  reached  the  plate  at  1.8  msec,  as 
computed  from  equation  2.2.  The  free  field  pressure 


transducer,  although  identical  to  the  other,  does  consistently 
give  earlier  signal  readings.  Table  V  indicates  that  a 
stand-off  distance  of  exactly  nine  feet  may  not  have  been 
maintained  between  the  charge  and  pressure  transducers. 

C.  POST-SHOT  EPSA  TESTING 

Initial  comparison  of  EPSA  predictions  and  experimental 
results  indicated  that  the  EPSA  model  was  inaccurate  or  had 
some  major  prediction  weaknesses.  However,  there  are  too 
many  variables  present  in  this  UNDEX  experiment  to  use  such 
a  gross  generalization.  The  mechanics  incorporated  in  EPSA 
are  sound  and  based  upon  singularly  proven  theories.  During 
an  experiment,  although  ideal  conditions  are  sought,  they 
are  seldom  totally  attained.  Adopting  a  more  realistic 
approach,  further  tests  were  conducted  using  EPSA.  The  main 
goal  of  the  post-EPSA  tests  were  to  try  to  achieve  the  measured 
deflections  using  the  code.  Computer  printouts  of  data  re¬ 
sults  and  PATRAN-G  color  slides  were  taken  of  the  various 
EPSA  runs,  but  they  are  too  lengthy  to  include  in  this  re¬ 
port.  Appendix  B  provides  sample  EPSA  input  decks  used  for 
the  different  runs,  should  there  exist  an  interest  in 
recreating  the  data. 

The  following  describes  a  list  of  post-EPSA  tests  that 
were  conducted  to  try  to  imitate  test  results.  First  a  9.64 
pound  TNT  run  was  made  but  the  deflections  were  not  signifi¬ 
cantly  different  from  that  of  an  eight  pounder.  Accounting 
for  the  fact  that  plate  tears  did  occur,  fundamental  fracture 


mechanics  were  attempted  to  be  incorporated  into  EPSA.  A 
change  to  the  subroutine  SM  STIF  allowed  for  modelling  the 
plate  with  the  stiffener  ends  cut.  This  provided  a  less 
rigid  test  panel  but  deflections  remained  symmetric.  The 
stiffener  ends  were  reattached  and  boundary  conditions  were 
written  into  the  EPSA  input  deck  to  simulate  actual  plate 
tears  along  the  top  and  side  edges.  The  results  from  these 
cuts  were  quite  interesting  as  this  was  the  first  time  non- 
symmetric  EPSA  deflection  patterns  were  produced.  However, 
these  cuts  did  not  portray  the  time  history  or  sequence 
whereby  these  tears  developed.  Finally,  using  data  points 
from  the  second  shot's  classic  free  field  pressure  plot, 
actual  pressure-time  history  information  was  inserted  into  a 
data  deck.  The  pressure-time  history  run  ignored  charge 
weight  and  tearing  conditions.  The  results  were  similar  to 
a  9.84  pound  charge.  This  was  expected  as  it  was  from  this 
graph  that  this  weight  was  calculated.  Using  Table  VI  and 
Figure  5.5,  a  complete  picture  of  EPSA  deflection  results 
are  given.  The  term  'cut'  in  Table  VI  refers  to  the  tear 
conditions  observed  from  the  second  shot.  None  of  the  above 
steps  reproduced  the  experimental  results,  but  they  did 
demonstrate  some  of  EPSA's  capabilities. 

D.  CONCLUSIONS  AND  RECOMMENDATIONS 

After  reviewing  the  results  of  this  experiment  and  the 
others  in  this  series,  it  is  clear  that  stiffener  tripping 
is  not  occurring.  Stiffener  tripping  is  not  experienced  due 


TABLE  VI 


Comparison  of  Actual  Deflection  to  EPSA  Predictions 


MARK 

ACTUAL 

8# 

9.84# 

0.16# 

P-T 

CUT 

8# 

CUT 

9.84# 

CUT 

P-T 

1 

.212 

.061 

.072 

.004 

.078 

.143 

.163 

.700 

2 

.583 

.328 

.366 

.031 

.406 

.618 

.681 

.679 

3 

1.19  5 

.619 

.676 

.063 

.688 

.958 

1.030 

1.030 

4 

.757 

.328 

.366 

.031 

.406 

_ 

5 

.211 

.061 

.072 

.004 

.078 

— 

_ 

lm 

6 

.448 

.138 

.156 

.017 

.115 

.294 

.328 

.329 

7 

.939 

.311 

.341 

.047 

.317 

.608 

.661 

.657 

8 

1.368 

.353 

.385 

.050 

.396 

1.150 

1.240 

1.220 

9 

1.042 

.311 

.341 

.047 

.317 

_ 

. 

10 

.524 

.138 

.156 

.017 

.115 

_ 

11 

.514 

.072 

.085 

.005 

.085 

.480 

.519 

.505 

12 

1.008 

.323 

.361 

.031 

.395 

.642 

.707 

.705 

13 

1.490 

.601 

.654 

.063 

.667 

1.010 

1.090 

1.090 

14 

1.082 

.323 

.361 

.031 

.395 

_ 

_ 

15 

.514 

.072 

.085 

.005 

.085 

_ 

_ 

16 

.474 

.061 

.072 

.004 

.078 

.438 

.484 

.473 

17 

.993 

.328 

.366 

.031 

.406 

.618 

.681 

.679 

18 

1.424 

.619 

.676 

.063 

.  688 

.958 

1.030 

1.030 

19 

1.025 

.328 

.366 

.031 

.406 

_ 

20 

.590 

.061 

.072 

.004 

.078 

_ 

21 

.314 

.138 

.156 

.017 

.115 

.133 

.152 

.159 

22 

.664 

.311 

.341 

.047 

.317 

.360 

.400 

.408 

23 

.944 

.353 

.385 

.050 

.396 

.515 

.556 

.  551 

24 

.713 

.311 

.341 

.047 

.317 

— 

25 

.383 

.138 

.156 

.017 

.115 

— 

_ 

26 

.106 

.099 

.114 

.011 

.088 

.097 

.112 

.122 

27 

.241 

.174 

.194 

.024 

.187 

.172 

.195 

.206 

28 

.334 

.172 

.194 

.024 

.194 

.236 

.256 

.258 

29 

.245 

.174 

.194 

.024 

.187 

_ 

30 

.129 

.099 

.114 

.011 

.088 

_ 

_ 

NOTE:  All  deflections  have  units  of  inches 

The  mark  number  refers  to  the  locations  in  Figure  5.5 


to  the  shearing  of  both  ends  of  the  stiffener.  Numerical 
analysis  does  not  predict  the  shearing  or  tearing  along  the 
sides.  From  the  data,  the  sequence  of  fractures  is  not 
discernible.  Knowing  the  order  of  fracture  occurrence 
would  give  a  better  clue  to  plate  responses.  Additionally, 
EPSA  predictions  do  not  accurately  forecast  experimental 
results.  In  general  the  code  predictions  are  lower  than 
measured . 

To  observe  stiffener  tripping,  the  stiffener  needs  to 
have  a  compression  mode  throughout  its  entire  length. 
Curvature  must  be  introduced  to  the  test  panel.  It  is 
recommended  for  future  studies  to  perform  UNDEX  testing 
using  a  stiffened  shell  model.  This  would  provide  all- 
around  compression  on  the  stiffener,  promoting  stiffener 
tripping  instead  of  shearing.  Also,  it  is  necessary  to 
incorporate  tearing  effects  or  fracture  mechanics  into  the 
numerical  analysis  scheme.  Symmetry  loading  conditions  are 
not  always  observed  from  an  UNDEX  experiment  and  these 
need  to  be  taken  into  account. 

It  is  recommended  that  UNDEX  testing  be  continued  at  NPS 
in  order  to  improve  experimental  techniques  and  to  acquire 
stiffener  tripping  data.  Once  tripping  has  been  accomplished 
then  investigations  can  branch  out  to  the  gross  shell 
responses  of  T  and  Z  stiffeners.  Improved  stiffener  tech¬ 
nology  would  lead  to  a  significant  advancement  in  ship 
design. 


2  3msec 


Field  Pressure  Gage--Shot  #1 


Figure  A. 4  Surface  Total  Pressure  Gage — Shot  #2 
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Figure  B.l  8  Pound  EPSA  Deck 
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Figure  B.2  Experimental  Pressure — Time  History 
Information 
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Figure  B.3  9.84  Pounds,  Cut — Top  and  Sides 
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